Numerical study of disorder effects on the three-dimensional Hubbard model 
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Combined effects of interactions and disorder are investi- 
gated using a finite temperature quantum Monte Carlo tech- 
nique for the three-dimensional Hubbard model with random 
potentials of a finite range. Temperature dependence of the 
charge compressibility shows that the Mott gap collapses be- 
yond a finite disorder strength. This is a quantum phase tran- 
sition from an incompressible phase to a compressible phase 
driven by disorder. We calculate the antiferromagnetic struc- 
ture factor in the presence of disorder as well. Strong anti- 
ferromagnetic correlation, which is characteristic of the Mott 
insulator, is destroyed by a finite amount of disorder. 
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I. INTRODUCTION 

As well as the repulsive interaction between electrons, 
the influence of disorder is essential in many electronic 
systems. Although either of these two effects causes 
a metal-insulator transition, the physical characters are 
quite different. The repulsive interaction tends to sup- 
press the double occupancy of electrons. On the other 
hand, in random potentials, electrons can favor doubly 
occupied states if the random potential is sufficiently low 
at the site. Therefore, the interaction and disorder may 
have opposite effects on the charge degree of freedom. 
In the insulator due to the interaction (Mott insulator), 
the charge fluctuation is strongly suppressed and a fi- 
nite charge gap opens, while the insulating phase due to 
disorder (Anderson insulator) does not necessarily have a 
charge gap. Another difference between these two insula- 
tors is the existence of magnetic correlation. Since the re- 
pulsive interaction induces local magnetic moments, the 
Mott insulator has a strong antiferromagnetic correla- 
tion. On the other hand, the Anderson insulator does not 
necessarily enhance magnetic correlation. Therefore one 
may expect that the interaction and disorder compete in 
both charge and spin degree of freedom. Especially, it is 
important to investigate to what extent the stability of 
the Mott insulator remains in the presence of disorder. 

The Hubbard model with disorder is one of the sim- 
plest models that include these two effects. In one and 
two dimensions, the transition from the Mott to the An- 
derson insulator is confirmed by various methods § - @. 
Also the dynamical mean field theory |l(J is applied to 
the infinite dimensional Hubbard model and consistent 
results are obtained jfl). On the other hand, there is 
no work beyond the mean field approximation in three- 
dimensional case WA. Therefore approximation- free re- 



sults can be useful to understand the three-dimensional 
strongly correlated system with disorder. 

In the present paper, we study the three-dimensional 
Hubbard model with random potentials using a finite 
temperature quantum Monte Carlo (QMC) method. The 
rest of this paper is organized as follows. In Sec. II, we in- 
troduce the three-dimensional disordered Hubbard model 
and describe physical observables. In Sec. Ill, we discuss 
the effects of disorder on the charge compressibility and 
the magnetic structure factor. 



II. MODEL AND METHOD 

The Hamiltonian of the disordered Hubbard model is 
given by 



(1) 



where t is the nearest neighbor hopping amplitude, 
is a nearest-neighbor link, U is the on-site interaction 
and {wi} are random potentials chosen from a flat dis- 
tribution in the interval [— W, W]. The system is on the 
cubic lattice in three dimensions and we use the periodic 
boundary condition. We treat the system in a grand 
canonical ensemble with the chemical potential fi. The 
grand canonical method is suitable to study the charge 
degree of freedom because the charge fluctuation is taken 
into account statistically |L3j. In this treatment, we make 
the system half-filled by setting /i = U/2. In the absence 
of disorder with sufficiently large U, the ground state is 
an antiferromagnetic insulator where the charge fluctua- 
tion is strongly suppressed. 

In order to obtain approximation-free results, we em- 
ploy a finite temperature determinantal Quantum Monte 
Carlo method Q, We also use the matrix- 

decomposition technique to remove numerical instabili- 
ties at low temperatures |lq ]. The simulations are per- 
formed in the half-filled sector (fi = U/2) for lattices with 
sizes uptoiV = 6x6x6 with U/t = 6. We choose a 
Trotter time slice size At ~ 0.15/t. We have checked 
that the systematic error due to the Suzuki- Trotter de- 
composition is almost independent of temperatures and 
does not change the qualitative feature. For each real- 
ization of disorder, we have typically run 2000 Monte 
Carlo sweeps for measurements after 500 sweeps in the 
warming up run. For all the observables, we average over 
24 realizations of disorder and the errors are estimated 
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by the variance among realizations of disorder. Since 
the system does not have a particle-hole symmetry in 
each realization of disorder, the negative-sign problem 
occurs. For example, the value of average sign is ~ 0.1 
for N = 4 x 4 x 4 and W/t = 1 at temperature T/t = 0.1. 
Although it is not so severe as a doped simula- 
tion at a very low temperature with strong disorder is 
difficult. 

The physical observables we have calculated are the 
compressibility k and the magnetic structure factor S(q) 
defined as 
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where iV is the number of sites, A e is the number of 
electrons and (3 is an inverse temperature. The charge 
compressibility k measures the charge fluctuation. If the 
system has a finite charge gap, the compressibility shows 
thermally-activated behavior in a low temperature region 
and vanishes at T = 0. On the other hand, the system 
without a charge gap has a finite compressibility at T = 
due to the existence of low-lying excitations. The mag- 
netic structure factor S(q) at q = (tt, tt, tt) diverges in a 
low temperature region when the system has an antifcr- 
romagnetic (quasi-) long-range order. 

III. RESULTS AND DISCUSSION 

Figure 1 shows the temperature dependence of the 
charge compressibility for several different strength of 
disorder. Average over 24 realizations of disorder is per- 
formed. Without disorder, the temperature dependence 
of the compressibility k shows thermally-activated be- 
havior reflecting the existence of a finite charge gap. The 
compressibility k at T = is zero within the numeri- 
cal accuracy for the pure system. It indicates that the 
ground state of the pure Hubbard model in three dimen- 
sion is in an incompressible phase. In the presence of 
disorder, the compressibility is enhanced. For weak dis- 
order, although enhanced, the compressibility still shows 
thermally-activated behavior and the value at T = 
seems to be zero. On the other hand, there is no ten- 
dency to decrease in k for W > W C (W C ~ U/2) down 
to the lowest temperature we studied. It means that the 
critical disorder strength to destroy the Mott gap, W c , 
is of the order of the Mott gap since the system is in 
the strong coupling region (U/t = 6). Although we can- 
not exclude the possibility of vanishing n at T = 0, what 
we have shown here is the best data within the numerical 
restriction. The results imply that sufficiently strong dis- 
order destroys the Mott gap which is of the order of the 
interaction in the strong coupling region. In the presence 



of disorder, when we discuss the physics locally, making 
one doubly occupied site gains an potential energy 2W 
at the maximum, while it costs a Coulomb energy U. 
Therefore one may expect that the Mott gap collapses 
at W > W C (W C ~ U/2) regardless of dimensionality in 
the strong coupling region. In other words, since charge 
properties of the Mott insulators in the strong coupling 
region is determined locally, the effect of disorder would 
be also local and independent of dimensionality. Indeed, 
the Mott gap collapses at W c ~ U/2 in one ||, two and 
infinite |l7| dimensions in the strong coupling region. It is 
in contrast to the quantum nature of the long-range prop- 
erties of the correlation functions which have a drastic 
difference in the dimensionality, (e.g. Luttinger liquid in 
one dimension) . The transition we observed is a disorder- 
driven quantum phase transition from an incompressible 
(gapped) to a compressible (gapless) phase. However, it 
does not necessarily mean an insulator-metal transition. 
The compressibility takes a finite value in both a metallic 
phase and an insulating phase due to disorder (Anderson 
insulator). It is possible that the competition between 
the interaction and disorder leads to a metallic phase es- 
pecially in the three-dimensional system. However, to 
distinguish these two phase, one needs simulations for a 
sufficiently large system, which we cannot perform be- 
cause of the negative-sign problem. 

Figure 2 shows the temperature dependence of the 
antifcrromagnetic structure factor S(tt,tt,tt). Since the 
ground state has an antiferromagnetic long-range order, 
the S(tt, tt, tt) shows diverging behavior toward the Neel 
temperature in the absence of disorder. For weak disor- 
der, the structure factor is slightly suppressed, but di- 
verging behavior is still observed down to the tempera- 
ture we studied. This means that the ground state still 
has an antiferromagnetic long-range order. When suffi- 
ciently strong disorder is included, the temperature de- 
pendence of the S(tt, tt, tt) changes qualitatively. The di- 
verging behavior of S(tt, tt, tt) disappears. This indicates 
that the long-range antiferromagnetic correlation is also 
destroyed by a finite amount of disorder. Ulmke et al. 
argue that weak disorder stabilizes antiferromagnetic or- 
der for U > U c , where U c is the interaction for which the 
Neel temperature takes a maximum value. , [fllf . Since 
the strength of the interactions we studied is U ~ U c |l8) , 
we did not observe it. 

In summary, the three-dimensional Hubbard model 
with random potential of a finite range has been studied 
numerically using a finite temperature quantum Monte 
Carlo method. The temperature dependence of the 
charge compressibility suggests that sufficiently strong 
disorder closes the Mott gap. The transition from an in- 
compressible phase to a compressible phase occurs at a 
finite strength of disorder. The disorder also destroys an 
antiferromagnetic long-range order which is characteris- 
tic of the Mott insulator. As in the case of the Mott gap, 
the antifcrromagnetic correlation is robust against weak 
disorder. These features are common in one, two and 
infinite dimensional systems. 
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FIG. 2. The antiferromagnetic structure factor as a func- 
tion of temperature (T/t), where U/t — 6, a)L = 4x4x4 
and b) L — 6 x 6 x 6. For weak disorder, the antiferromag- 
netic structure factor shows diverging behavior down to the 
temperature we studied. On the other hand, for strong disor- 
der, the divergence behavior is not observed. This indicates 
that sufficiently strong disorder destroys the long-range anti- 
ferromagnetic correlation which is characteristic of the Mott 
insulator. 
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FIG. 1. Temperature dependence of the charge compress- 
ibility k, where U/t = 6, a)L = 4x4x4 and b) L = 6 x 6 x 6. 
Without disorder, k shows thermally-activated behavior and 
decreases toward T = 0, indicating the existence of a charge 
gap. For weak disorder, « still shows thermally-activated be- 
havior. On the other hand, for strong disorder, « does not 
decrease down to the temperature we studied. It is a disor- 
der-driven transition from an incompressible phase to a com- 
pressible phase. 
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